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Abstract: 

The geometric clustering of partons in the transverse plane of nuclear collisions leads 
^ . for increasing A or y/s to percolation. In the resulting condensate, the partons are de- 
confined but not yet in thermal equilibrium. We discuss quarkonium dissociation in this 
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precursor of the quark-gluon plasma, with an onset of dissociation when the saturation 
scale of the parton condensate reaches that of the given quarkonium state. 
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O , Statistical QCD predicts color deconfinement for sufficiently hot strongly interacting 

^ I systems in full equilibrium. In the resulting quark-gluon plasma, both the momenta 

■ and the relative abundances of quarks, anti-quarks and gluons are determined by the 

p ■ temperature of the medium. It is not evident if and at what evolution stage high energy 
^ . nuclear collisions produce such equilibrium, nor is it evident that color deconfinement is 
D ! restricted to such ideal thermal systems. It thus seems natural to ask what conditions are 
^ ' necessary in the pre-equilibrium stage to achieve deconfinement and perhaps subsequent 
quark-gluon plasma formation. In recent years, the occurrence of color deconfinement in 
^ ■ nuclear collisions without assuming prior equilibration has therefore been addressed on 
the basis of two closely related concepts, parton percolation |^ and parton saturation 

iii- 

Both start from the observation that in a central nucleus- nucleus collision at high 
energy, one finds in the transverse nuclear plane interacting parton^ of different transverse 
scales. At low densities, one can define individual partons originating from nucleons of 
the incident nuclei. Once the density of partons becomes so high that they form a dense 
interacting cluster, independent parton existences and origins are no longer meaningful: 
the resulting cluster forms a condensate of deconfined partons. The condensate is formed 
in the sense of droplets condensing to form a liquid, and the partons which make up this 
condensate are no longer constrained by any hadronic conditions. 

Consider a distribution of partons of transverse size vrr^ over the transverse nuclear 

relation between deconfinement and percolation was suggested quite long ago Q; but the first work 
on color percolation in nuclear collisions was given in terms of strings rather than partons. We shall 
nevertheless restrict ourselves here to a parton picture [0 . 



plane nR^, with r << The fundamental aim of parton percolation studies [0, 0, ^ is 
the determination of the transition from a normal hadronic collision situation of disjoint 
partonic 'discs' to a connected cluster of such discs spanning the nucleus, the parton 
condensate. Percolation theory predicts this transition to occur when 

— ^ 7rr^ = Uq Tcr"^ = t] ^ ijc = 1.128 (1) 

i.e., when the parton density r] = n^Trr^ measured in terms of parton size reaches the 
percolation threshold rjc. In the 'thermodynamic' limit of infinite spatial size Ra — ^ oo 
and infinite parton number Nq —>■ oo, the largest connected cluster first spans the system 
at this point. The critical value rjc = 1.128 is determined by extensive numerical studies 
0. Since the partons overlap, this does not mean that the entire transverse nuclear 
surface is covered by parton discs. In fact, at the percolation point, only the fraction 
1 — exp(— r/c) — 2/3 of the nuclear area is covered by partons. At rj = rjc, the critical 
clustering behavior of the system can be specified in the usual way in terms of critical 
exponents. In particular, the size S{ri) of the largest cluster diverges for r] ^ rjc as 

S{v)-^{Vc-v)-\ (2) 

with the critical exponent 7 = 43/18. While this holds strictly only for infinite systems, 
it is verified that even for rather small spatial systems, the transition from very small 
size to percolating cluster occurs in a very narrow density interval. In other words, even 
at finite size there is almost critical behavior. This will become quite important for our 
further considerations. 

The essential idea of parton saturation is that the increase in the number of partons 
for small x, as obtained from deep inelastic scattering experiments, must stop when the 
density of partons becomes so high that they overlap and form large interacting clusters; 
fusion and splitting then causes their number to approach a constant. The onset of 
saturation has been discussed in various ways; making use of their transverse size, it can 
also be quite naturally determined by the percolation condition (|l]), which then fixes the 
saturation scale in terms of A and the c.m.s energy ^/s. Saturation sets in when the 
parton density Uq in terms of the partonic interaction cross section aq approaches the 
critical value rjc = 1.128. The partonic cross section depends on its inverse transverse 
momentum fc-r, (^gikx) ~ I/^t? ^^"^ ^^e parton density ng{k^) for a fixed resolution scale 

can be obtained from the gluon distribution function determined in deep inelastic 
scattering. The novel aspect from the point of view of percolation is that the density of 
partons nq{k'^) is related to their transverse size a{k'^), so that with the functional form of 
these two quantities given, the percolation condition specifies the scale of the percolating 
partons. Let us consider this in more detail. 

The distribution of partons in an incident nucleon of momentum p is given in terms of 
their fractional momentum x c::^ k/p and the resolution scale Q^, through integration over 
the partonic transverse momentum kr- The relevant resolution scale in a nucleon-nucleon 
collisions is the largest transverse momentum fc™*^^ = Q for which partons are resolved. 

^For simplicity, we assume the partons to have a fixed transverse radius; the extension to a distribution 
of radii is straight-forward and does not change the picture. 
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The number of partons at fixed x and integrated over kx < Q is given by the sum of the 
contributions of gluons plus those of sea quarks and antiquarks, 



dy 



xg{x, Q"^) + Y^[xqi{x, Q^) + xqi{x, Q^)], 



(3) 



where g{x,Q'^) denotes the gluon distribution function, qi{x,Q'^) and qi{x,Q'^) that of 
up and down (z = 1,2) quarks and antiquarks, respectively. The distribution functions 
are determined from parametrisations of deep inelastic scattering data, and thus eq. (H) 
provides the number of partons at central rapidity y = 0, with x = Q/y/s. We shall see 
shortly that for RHIC and higher energies, the gluon contribution is strongly dominant; 
for SPS energy, however, the quark and antiquark contributions cannot be neglected. 

We further need the parton size. In the simplest percolation approach, which we shall 
follow here, this is just the geometric cross section vr/Q^. More dynamical considerations 
lead to numerical modifications, ag{Q^) = Kas{Q'^)'n/Q'^, where as(Q^) is the running 
couphng at scale and k a given constant. With the geometric cross section we obtain 
that in an A A collision at y = 0, the equation 



determines the onset of percolation. Here ns{A) specifies the density of parton sources 
in the transverse plane of a central AA collision. At SPS energy, this is essentially the 
density of wounded nucleons [§], ns{A) ~ n^{A). For higher energies and harder partons, 
collision-dependent contributions will play a significant role, and so a more suitable form 
here is a combination of the two sources. 



with < X < 1 0. 

As determined by eq. (^, partons with kx < Qs condense to form an overlapping 
and hence interacting cluster spanning the system, the parton condensate. Within this 
cluster, they can fuse or split and thus lose their independent existence. We recall that 
at this point, 2/3 of the nuclear area is covered by the parton condensate. 

The relation between saturation and percolation has so far not been much empha- 
sized. For a study of the new percolating medium at very large A or y/s, the 'color glass 
condensate' of Ref. it is indeed not so important. However, it does become crucial 
for a detailed picture of the onset of parton condensation. We know from percolation 
theory that in the large volume limit this is a critical phenomenon and hence even for 
finite systems takes place in an almost singular way. 

To study the onset of parton percolation in AA collisions, it is convenient to rewrite 
eq. @ in the form 




ns{A) = {1- x)n^ + XTicoii, 



(5) 




(6) 
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with the density of parton sources Ug in fm~^; the factor 25 arises when this is converted 
to GeV^. In eq. (|), the hadronic parton distribution {dNq/ dy) /Q"^ is compared to the 
density of parton sources in a nucleus-nucleus collision. For low source densities, i.e., for 
small A, 8.98/ns(y4) remains well above the A-independent l.h.s., which is a function only 
of Q and ^/s^. there are not enough partons to form a condensate. For sufficiently large 
A, however, 8.98/ns(v4) intersects {dN^/ dy) /Q"^ at some Qs, thus defining the onset of 
percolation. All partons with < Qg merge to form the condensate, in which interactions 
prevent much further increase in the number of partons, i.e., there is saturation. 

To see when that happens, we have to make use of a specific set of parton distribution 
functions. The kinematic range relevant for our analysis is 0.5 GeV^ < < 2 GeV^, 
with 0.02 < a; < 0.1 for SPS and 0.003 < a; < 0.02 for RHIC. Among the commonly used 
PDF parametrizations, only the set GRV94 [|1^ goes down to such small values of Q^. 



This parametrization describes well the available data on the proton structure function 
from the E665 and NMC collaborations for 0.4 GeV^ < Q"^ < 2 GeV^ and x < 0.01 
111 , |12|. It also reproduces quite well the small x HERA data |]10[ |13] for x > 10~^. We 



have therefore calculated (dNg/dy) using the next-to-leading order GRV94 PDF's in the 
DIS scheme (GRV94DI) [|10|; the resulting (dNg/dy) /Q'^ are shown in Figs, la and lb for 
SPS (20 GeV) and RHIC (200 GeV) energies, respectively. 

Any uncertainty in (dNg/dy) is mainly due to the gluon distribution. At SPS energy, 
this can be estimated by comparing calculations using leading and next-to-leading order 
GRV94 PDF's. The resulting uncertainty is below 3 % for > 0.6 GeV^ and increases 
to at most 10 % at 0.5 GeV^. At RHIC, the next-to-leading order GRV94 PDF's are 
consistent with the gluon distributions determined by the HI and ZEUS collaborations 
and with the constraints from charm measurements fl^, while the leading order results 
for the gluon distribution are too large. This gives us some confidence in using the next- 
to-leading order GRV94 parametrization. 

We also note that these gluon distributions are not very different from those proposed 



in more recent phenomenological saturation studies for x < 0.01 [15|. This approach is 
very successful in describing the low x HERA data and it can in fact also account for the 
E665 data on F2 in the kinematical region relevant for RHIC. The gluon distribution of 
jTSf is directly related to Ff; the sea quarks are present only virtually as small dipoles from 



photon splitting |]I6[. In such an approach, the uncertainty in the gluon distribution can 
be avoided, and it has been used to predict the A/s-dependence of hadron multiplicities 
at RHIC 0. 

Once the density of the percolating medium is sufficiently high, the PDF approach of 



Ref. to the gluon distribution is more appropriate to study the resulting condensate; 
however, it is not suitable to study the onset of condensation. In this connection we note 
that the effect of the sea quarks cannot be completely neglected even at RHIC energy. 
In general, the relative weight of the different parton species in the produced condensate 
will depend on Q and a/s, and with increasing energy, the condensate becomes more and 
more gluon-dominated. Thus the ratio of gluons to quarks and antiquarks in the parton 
condensate based on the GRV94 PDF's is at Q = 1 GeV found to be about 1.2 for SPS 
energy and 4.0 for RHIC; in a chemically equilibrated quark-gluon plasma, it is about 0.5. 



Hence before any thermalization, the medium is strongly gluon-dominated |19 . 
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In order to illustrate the effect for different central AA collisions, we assume a spherical 
nuclear profile, which gives at SPS energy, with x = in eq. (H), 

2A 

n,{A)=n^{A) = -^. (7) 

Inserting this into eq. (|^), together with the values of [dNq{x,Q'^)/dy]/Q^ at y/s = 20 
GeV of the previous section, we obtain the results shown in Fig. la. The intersection 
points determine the onset of percolation. It is seen that in this very simplified picture, 
parton condensation begins for A ~ 60, with Q ~ 0.7 GeV. To obtain the corresponding 
behavior at -y/i = 200 GeV, nucleon collisions have to be included as source of partons. 
With 

3 / A*^/^ \ 

and X = 0.09 in eq. (|), we find the percolation points in Fig. lb. In eq. (|]), the factor 
3/4 comes from averaging over the nuclear profile. In Fig. 2, the percolation values Qg 
for central AA collisions are displayed as function of A. At SPS energy, we thus do not 
obtain parton condensation below A ~ 60; at RHIC energy, the higher parton density 
lowers the onset to A ~ 40. 

In the experimental study of J/ip suppression, the production is measured at different 
centralities, so that we now have to determine the parton source density at fixed A and 
varying impact parameter. This is done in a Glauber analysis based on Woods-Saxon 
nuclear profiles, with a collision-determined weight [l^ . In Fig. 3a, we show the resulting 



percolation behavior as function of the effective number A^part of participants in a Pb — Ph 
collision at ^/s = 20 GeV. The threshold for parton percolation is found to be slightly 
below A^part — 150. The corresponding calculations for Au — Au collisions at ^/s = 200 
GeV, with a collision-dependent term in and again x = 0.09 in eq. (|^), lead to the 
results shown in Fig. 3b. The onset of parton condensation at RHIC is thus shifted 
to considerably more peripheral collisions. In Fig. 4, the centrality dependence of the 
percolation scale Qs is shown; at the onset point, the condensate contains partons of 
different sizes r, with r > 1/(0.7 GeV) at SPS and r > 1/(0.9 GeV) at RHIC. 

We now turn to the effect of parton condensation on J/ip production. It is known 
from pA collisions that normal nuclear matter leads to reduced charmonium production. 
Therefore we first have to consider such 'normal' suppression, since for pA collisions, at 
least up to RHIC energies, we are well below the threshold for parton condensation^. 
Consider the production of charmonium in the nuclear target rest frame. A gluon from 
the incident proton fluctuates into a virtual cc pair; in a collision with one of the target 
nucleons this is brought on-shell, its color is neutralized, and it eventually becomes a 
physical J/ip of size dj/^ ~ 0.4 fm. Thus in the rest-frame of the cc pair, a time of at 
least 0.4 fm is needed for J/ip formation. During its pre-resonance stage, for t^c < 0.4 fm, 
it will travel a distance 

SXp 



d - = T 



4M2 

1 + Wl + ^ 



(9) 



determination of the collision energy for which parton saturation occurs in ppjpp or pA collisions 
is presently difficult, since it involves .T-values much smaller than presently accessible in deep inelastic 
scattering. Disregarding this small x-behavior leads to a very early onset of parton percolation p2[. 
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in the target rest frame; here m denotes the mass of the nucleon, M that of the cc pair 
and xp the cc Feynman momentum fraction. For the production of a cc at rest in the 
nucleon-nucleon c.m.s., this becomes dec = Tcc{\^/2m). From this it is immediately seen 
that at ^/s ~ 40, and xp > 0, the nascent J/ip effectively traverses even heavy nuclear 
targets in its pre-resonance stage. The situation is very similar for Xc and ip' production, 
for which the pre-resonance life-times are if anything even larger. Thus for the mentioned 
^/s and xp, the nuclear target sees of all charmonium states only the small pre-resonance 



precursor, so that all should suffer the same degree of nuclear suppression ||2^. For 
negative xp (and perhaps to some extent also for lower a/s), the nucleus should begin 
to see physical resonances, and as a result the suppression should become larger for the 



higher excited states with their larger radii ||2T 



On the basis of this information, normal charmonium suppression has generally been 
studied in terms of pre-resonance dissociation in standard nuclear matter, leading to a 
break-up cross section around 5-6 mb This implies a mean free path of about 



12 fm. For a multiple collision analysis of Glauber type, the mean free path has to 
exceed the coherence length (essentially the size) of the pre-resonance in the rest frame 
of the nucleus: otherwise, the cc is wounded several times before it has had a chance 
to register the first interaction. For collision energies ^/s < 40 GeV, i.e., in the range 
of fixed target experiments, this condition is satisfied, with coherence lengths below 8 
fm. At RHIC energy, on the one hand the cross section for pre-resonance break-up in 
normal nuclear matter could change, and on the other hand the coherence length becomes 
dilated ten times more. Now interference effects of the Landau-Pomeranchuk-Migdal type 
may have to be taken into account, which could lead to a reduction of 'normal' nuclear 
suppression |2^. Hence measurements of charmonium production in pA collisions at 



RHIC are absolutely essential for an understanding of whatever J/ip suppression in AA 
collisions is observed there. 

We now want to consider the additional 'anomalous' suppression of charmonium pro- 
duction due to the dense partonic medium created in AA collisions. The virtual partons 
in the incoming nuclei coalesce to form a condensate in a time tc — 1/Qs determined by 
the saturation scale; in the color glass picture 0, this is the time needed to melt the 
frozen glass. The interacting and expanding parton condensate can subsequently lead to 
the formation of a thermalized quark-gluon plasma; a crucial factor for this is the energy 
density reached at thermalization. We restrict ourselves here to the parton condensate 
stage. 

In AA collisions at RHIC energy, the colliding nuclei are in the overall c.m.s Lorentz- 
contracted to about 0.1 fm. They will therefore sweep past the nascent charmonium in 
its pre-resonance state and before parton condensation sets in, resulting in some form 
of normal nuclear absorption. After about 0.2 fm, the nuclei are well out of the way 
and the parton condensate is formed. Any produced and surviving charmonium states 
now encounter this new medium, either as fully formed resonances or in the late pre- 
resonance stage. For SPS energy, a similar discussion is more complex, since the nuclei are 
contracted to only about 1 fm diameter; this will introduce a smearing in the comparison 
of the different time scales. We shall here neglect this and assume that also at the SPS 
there is first pre-resonance absorption in normal nuclear matter, followed by the effect of 
the parton condensate on the survivors. 
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In the first attempt to describe J/ip suppression in terms of parton percolation 0, 
it was assumed that different charmonium states i define particular scales Qj, and the 
onset of percolation for partons of that scale leads to the dissociation of all charmonia 
of that species within the percolating cluster. The number of partons was taken as 
a scale-independent function of ^/s. Parton saturation in fact specifies [dNq/ dy\{Ql) as 
function of the scale Qs-, so that we now have (see Figs. 2 and 4) at given collision energy a 
parameter- free determination of the onset line Qs(A) or Qs{Npart) of parton condensation. 

Within the parton condensate, color fields of a Q^-dependent strength will affect the 
binding of a cc dipole charmonium state. This effect can be addressed in different ways. 
One possible way is to consider the cc propagation in a classical field, which can represent 
the gluon field of the condensate or that of the nucleus; the information on specific 
medium is encoded in the corresponding field correlator [^. For a small singlet dipole 
the probability to remain singlet is proportional to exp(— /tr^), where r is the dipole size 
and K is a dimensional parameter determined by the field strength. In the color glass 
approach, k ~ Thus for small singlet dipoles, Qlr"^ « 1, the probability to 

remain in a singlet state is close to unity (color transparency). Motivated by this fact 
we shall here adopt the simple model of Ref. 0, assuming that a charmonium state i 
of scale Qi will be dissociated if it finds itself in a parton condensate of scale Qs ^ Qf, 
otherwise it will survive. This very simplistic picture allows an analysis of nuclear profile 
effects (condensed and non-condensed regions of the collision profile) and thus provides 
some direct predictions for the centrality-dependence of anomalous J/ip suppression for 
given A and -y/i. Other approaches that have been suggested include the study of the 



time evolution of the screening masses in the parton cascade |25] and in the color glass 



condensate 26 



The radii of the observable charmonium states as obtained from the solution of the 



Schrodinger equation with Cornell potential are 27 



rj/^ ^ (0.9 GeV)-\ ~ (0.6 GeV)-\ - (0.45 GeV)~^ (10) 

These have to be compared to the partonic saturation radii shown in Fig. 4 as 
functions of the number of participants in central AA collisions. It is seen that at SPS 
energy, the onset of Xc and ip' suppression effectively coincides with the onset of parton 
condensation, while the J/ ip survives up to larger Qs and higher parton densities^. Hence 
in Pb—Pb collisions, anomalous suppression should start with the elimination of J/ip feed- 
down from Xc and ip' at Npart — 150; the further suppression of direct J/ip production 
should set in for Npart — 250. In Fig. 5 we show the suppression pattern observed by the 



NA50 collaboration at the CERN-SPS [28|; the two 'steps' observed in the anomalous 



suppression pattern agree fairly well with the expected onsets of parton condensation. 
Note that at the percolation points, the percolating cluster covers only a fraction of the 
transverse area; with increasing centrality, this fraction increases. Hence more detailed 
studies are needed to determine the actual amount of suppression as function of centrality 
(see 0); the present study only indicates the onset points. 

*The fate of the ip' is most likely more complex. It should certainly be dissociated once parton 
condensation sets in; however, it is very much more weakly bound than Xc smd J/ijj, so that a less dense 
environment could also lead to its break-up. 
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Finally we want to consider the suppression pattern expected for RHIC experiments. 
From Fig. 4, it is seen that for Npart > 80, all charmonium states should suffer anomalous 
suppression, so that here there should only be one onset point. Moreover, collisions 
with Npart = 80 correspond to an impact parameter of 6 ~ 10 fm, and this may be 
too peripheral for meaningful measurements. For a study of the onset of anomalous 
suppression at RHIC, it will thus most likely be necessary to study AA collisions for 
much lower A, as already noted previously [Q. 

In conclusion, we have shown that a very simplistic parton percolation approach leads 
to a conceptually reasonable and effectively parameter-free description of the observed 
anomalous J/ip suppression pattern. 
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Figure 3: Parton distribution vs. with the percolation hmits for different centralities, 
given by the number of participant nucleons; : (a) Pb — Pb at y/s — 20 GeV, (b) Au — Au 
at = 200 GeV. 
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Figure 4: Saturation momenta Qs vs. AA for central Ph — Pb collisions at ^/s = 20 GeV 
and central Au — Au collisions at \/s = 200 GeV, together with the critical scales for 
anomalous J/i^, Xc and ip' suppression. 





1.4 r 




1.2 - 




1 - 


o 


0.8 - 


(/) 




Cl 




X 
CD 


0.6 - 


CO 






0.4 - 




0.2 - 



-//- 



>K p-p, p-D 

• p-Cu, p-W, p-U (200 GeV) 

□ p-C, p-AI, p-Cu, p-W (450 GeV) 

-//- 



0-Cu, 0-U 
S-U 

Pb-Pb 1996 
Pb-Pb 1996 (min.bias) 
Pb-Pb 1998 (min.bias) 



i. 



<; 



2 3 4 5 



50 100 150 200 250 300 350 400 
Number of participants 



Figure 5: The J/ip survival probability determined at SPS energy as function of the 
number of participant nucleons, from ||28| . 
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